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We report on the infrared limit of the quenched lattice Landau gauge gluon and ghost propagators 
computed from large asymmetric lattices. In particular, the compatibility of the pure power law 
infrared solutions of the Dyson-Schwinger equations with the lattice data is investigated and the 
exponent k is measured. The gluon lattice data favour K ~ 0.52, which would imply a vanishing 
zero momentum gluon propagator. For the subset of lattices where the ghost propagator was 
computed, the data are not compatible with a pure power law. Our data also show a decreasing 
running coupling in the infrared region. Furthermore, positivity violation for the gluon propagator 
is also verified. 
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1. Introduction and motivation 



The study of the infrared limit of Quantum Chromodynamics (QCD) requires the use of non- 
perturbative methods. Two first principles non-perturbative approaches are Dyson-Schwinger equa- 
tions (DSE) and the lattice formulation of QCD. Because both methods have good and bad features, 
a comparison between their results is a good test of our understanding of the low energy limit of 
QCD. During the last years, there has been an effort to compute the infrared gluon and ghost prop- 
agators in Landau gauge, 



<(*) = S^(V-^)»(<7 2 ), (1.D 

G ah {q) = -8 ah G{q 2 ), (1.2) 

and a running coupling constant 1 defined from these propagators 

<Xs{q 2 ) = a s {n 2 )Z 2 ghost (q 2 )Z gh<on (q 2 ); (1.3) 

Z g host{q 2 ) = q 2 G(q 2 ) and Z g i uon {q 2 ) = q 2 D(q 2 ) are the ghost and gluon dressing functions. 

In [[I]], assuming ghost dominance, it was computed a solution of the DSE that predicts pure 
power laws for the propagators, namely 



Z gluo „(q 2 ) ~ (q 2 ) 2K ,Z ghos! (q 2 ) ~ (q 2 )- K , (1.4) 

with K = 0.595. This implies a vanishing (infinite) gluon (ghost) propagator for zero momentum. 
Other studies of the infrared limit also predict K > 0.5 @ ^[ §]. 

As an infrared analytical solution of the DSE, the pure power laws are valid only for very low 
momenta. Indeed, comparing the DSE solution for the gluon propagator [||] with the corresponding 
pure power law, see figure p], it comes out that the power law is valid only for momenta below 200 
MeV. 

For lattice QCD, it is a challenge to perform a simulation with a minimum number of points 
in the region of interest. The symmetric lattices available at the moment have a limited number of 



points in the infrared region (see, for example, [ 13 , |22|, |25j] ) 



In the series of papers ^ [7| [| K], 1 1] we have been using large asymmetric lattices xL t , 
with L t S> L s , to investigate the infrared limit of the gluon and ghost propagators. In this article 
we report on the status of our investigations concerning the use of asymmetric lattices to study the 
infrared properties of QCD. 



2. Gluon propagator 

In [0], we have computed the gluon propagator for SU(3) four-dimensional asymmetric lattices 
L A 3 x256, with L s = 8, 10,..., 18. 



See [1230 for a recent review on the running coupling constant. 
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Figure 1: The gluon DSE solution [||] compared with the corresponding pure power law. 
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Figure 2: On the left, the gluon propagator for 16 3 x 128 and 16 3 x 256 lattices, considering only pure 
temporal momenta. Note the logarithmic scale in the vertical axis. On the right, the gluon propagator for all 
lattices L 3 x 256. For comparisation, we also show the 16 3 x 48 and 32 3 x 64 propagators computed in 



As reported in [||, ||, [R], [TT, 12, [13|], there are clear finite volume effects. However, the various 
simulations performed up to now show that the approach to the infinite volume limit L s — > +oo is 
smooth. 

In what concerns the finite volume effects, our simulations suggest that the finite volume 
effects are essentially due to a relatively small spatial extension. Indeed, the gluon propagator 
was computed for 16 3 x 128 and 16 3 x 256, and the data are undistinguishable — see figure || 
(left). This result gives us confidence that the temporal size of our lattices is sufficiently large. 

However, it was observed that the propagator depends on the spatial size of the lattice — see 
figure |2| (right). The gluon propagator decreases with the volume for the smallest momenta and 
increases with the volume for higher momenta. 



2.1 Infrared exponent 

In order to compute the infrared exponent K from the lattice data, we considered fits of the 
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smallest temporal momenta of the gluon dressing function Z g i uon (q 2 ) = q 2 D(q 2 ) to a pure power 
law with and without polinomial corrections 2 . In general, the K values increase with the volume of 
the lattice. So, our K can be read as lower bounds in the infinite volume figure Koo. 

The gluon propagator can be extrapolated to L s — > +00, as a function of the inverse of the vol- 
ume, fitting each timelike momentum propagator separately, and assuming a sufficient number of 
points in the temporal direction. Several types of polinomial extrapolations were tried, using differ- 
ent sets of lattices, and we conclude that the data are better described by quadratic extrapolations 
of the data from the 4th and 5th largest lattices. 

The values of k extracted from these extrapolated propagators are K = 0.5215(29), with a 
X 2 jd.o.f . = 0.02, using the largest 5 lattices in the extrapolation, and K = 0.4979(66), % 2 /d.o.f. = 
0.27 using the largest 4 lattices. Fitting the extrapolated data to the polinomial corrections to the 
pure power law, one gets higher values for K. The first value K = 0.5215(29) is on the top of 
the value obtained from extrapolating directly fc as a function of the volume [gj. Note also that in 
[^, [TTJ], we fitted the gluon temporal data for larger ranges of momenta, using other model functions, 
giving always values for K above 0.5, supporting again an infrared vanishing gluon propagator. In 
conclusion, one can claim a K £ [0.49,0.53], with the lattice data favouring the right hand side of 
the interval. 



2.2 Positivity violation 

In QCD, the violation of reflection positivity for the Landau gauge gluon propagator means 
that the gluon cannot appear as a free asymptotic S-matrix state. This may be viewed as an indica- 
tion of gluon confinement. 

On the lattice, one can study positivity violation for the gluon propagator from the real space 
propagator, 

/oo 
dpD(p,p = 0)exp(-ipt). (2.1) 
-00 

Finding C(t) < for some t is a sign of positivity violation. 

Positivity violation for the gluon propagator has been observed for the Dyson-Schwinger gluon 
propagator [14], as well as for the propagator computed from symmetric lattices [24, |13|]. Here we 
show the real space propagator computed from our asymmetric lattices, including in figure [|, C(t) 
for the two infinite volume extrapolations of the gluon propagator computed in [§]. 

We thereby confirm that positivity violation occurs for the gluon propagator. Our data show 
that the time for positivity violation to happen, decreases when the spatial lattice volume is in- 
creased. Furthermore, the infinite volume limit suggests that positivity violation shows up at 
t ~ \.5fm. Previous studies show similar values [13, [l4|]. Moreover, similarly to what was ob- 
served in [24], for large time separations our data show an oscilatory behaviour. 



3. Ghost propagator 



So far, we have computed the ghost propagator for our smallest lattices [10] (10 3 x 256, 12 3 x 
256, 16 3 x 128), using both a point source method [ 18] (we averaged over 7 different point sources 



2 The results can be seen in table II of [ 
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Figure 3: On the left, the real space gluon propagator C(t) for our largest lattices and for the two extrapola- 
tions considered in [Q]. On the right, a zoom of the region of interest showing all lattices. 



to get a better statistics), and a plane-wave source [15]. In both cases we used the pre-conditioned 
conjugate gradient algorithm, as described in [|lq]. The plane- wave source method provides better 
statistical accuracy, but we can only obtain one momentum component at a time. The point source 
method allows to get all the momenta in one go, but with larger statistical errors. 

In figure |] (left) it is shown the ghost dressing function Z^ wst (q 2 ) = q 2 G(q 2 ) for the 16 3 x 128 
lattice, gauge fixed using CEASD method [O], computed with both methods. As in the gluon 
case, we can see differences, in the infrared, between pure temporal and pure spatial data. These 
differences vanish for sufficiently high momenta. 
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Figure 4: On the left, the bare ghost dressing function for the 16 3 x 128 lattice, gauge fixed with CEASD 
method. "p2p" ("pws") stands for the ghost components computed using a point (plane wave) source. On 
the right, the bare ghost dressing function computed from a plane wave source for all lattices. 



On the right hand side of figure one can see the ghost dressing function only for the plane- 



wave data, for the available lattices. As in the gluon case []11[], we are able to evaluate the effect 
of Gribov copies on the lattice 16 3 x 128 by considering different gauge fixing methods. The data 
show clear effects of Gribov copies over a large range of momenta, as expected from other studies 
[£5],[l6|]. Also, we see finite volume effects if one compares propagators from lattices with different 
spatial sizes. 

In what concerns the infrared region, we were unable to fit a pure power law, even considering 
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Figure 5: On the left, the running coupling constant for the 16 3 x 128 lattice, gauge fixed with CEASD 
method. On the right, the running coupling constant for all lattices. The data were computed using a plane 
wave method. 



polinomial corrections [10]. 



4. Running coupling constant 

From the gluon and ghost dressing functions, one can define a running coupling constant 
see eq. 



1.3 



The DSE infrared analysis predicts a running coupling constant at zero momentum different 
from zero, aj(0) = 2.972 [p. On the other hand, the DSE solution on a torus [19], and results from 

show a decreasing coupling constant for small momenta. Using an 



lattice simulations [ |16| , [21[ |1( 
asymmetric lattice allow us to study smaller momenta having in mind to provide, at least, a hint to 
this puzzle. 

Again, our lattice data show finite volume effects, if one compares pure temporal and pure 
spatial momenta, see figure |5| (left). Comparing the results for all available lattices (plane-wave 
source), see figure || (right), we can see, as in the ghost case, finite volume effects, and clear Gribov 
copies effects. 

In what concerns the infrared behaviour, we tried to fit the lowest momenta to a pure power 
law, {q 2 ) Ka . We concluded that this power law is only compatible with the data from 16 3 x 128 
lattice, gauge fixed with CEASD method, giving K a ~ 0.688, with % 2 /d.o.f. ~ 0.011. The reader 
should be aware that it is also possible, in some cases, to fit the infrared data to a (0)(1 + aq 2 + . ..) 
and get a a(0) / 0. Therefore, we can not give a definitive answer about the behaviour of the 
running coupling constant for q = 0. Note, however, that OCs(q 2 ) for the smallest momenta, seems 
to increase as a function of the volume. 



5. Future work 



Currently, we are engaged in improving the statistics for our larger lattices and the extrapo- 
lations to the infinite volume limit for the gluon and ghost propagators. Furthermore, we plan to 
perform simulations with larger lattices and, hopefully, combine all our results to provide a reliable 
answer on the behaviour of the infrared QCD Green's functions. 
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